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Abstract 

We demonstrate the feasibility of a nonperturbative analysis of quantum field 
theory in the worldline formalism with the help of an efficient numerical algorithm. 
In particular, we compute the effective action for a super-renormalizable field theory 
with cubic scalar interaction in four dimensions in quenched approximation (small- 
iVf expansion) to all orders in the coupling. We observe that nonperturbative effects 
exert a strong influence on the infrared behavior, rendering the massless limit well 
defined in contrast to the perturbative expectation. Our numerical method is based 
on a direct use of probability distributions for worldline ensembles, preserves all 
Euclidean spacetime symmetries, and thus represents a new nonperturbative tool for 
an investigation of continuum quantum field theory. 

PACS numbers: ll.10.Gh, ll.I5.Tk, I2.20.Ds 



1 Introduction 

Many current problems in particle physics as well as statistical physics demand for new 
nonperturbative field theoretic techniques. Particularly, the treatment of strongly coupled 
fluctuations requires methods that go beyond standard perturbative techniques. 
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Conventionally, perturbation theory is understood as an expansion for both small cou- 
pling and small amplitude of the fields. The case of larger fields but small coupling can 
be addressed with effective-action techniques that take the coupling to all orders to an 
external field into account; prominent examples are the Heisenberg-Euler effective action 

or the Coleman- Weinberg effective potential [2]. 

In the present work, we propose a technique that has the potential to deal with the 
opposite limit: arbitrary coupling but weak external-field amplitude 1 . The technique is 
based on the worldline formulation of quantum field theory that goes back to ideas of 
Feynman it can be viewed as a mapping of field theoretic problems onto quantum 
mechanical path integrals, identifying the trajectory of a fluctuation in coordinate space 
with the path of a quantum mechanical particle. The approach had occurred rarely in the 
literature, see, e.g., j3], until it was realized that certain perturbative computations simplify 
tremendously in the worldline approach E] ; the deep reason behind this observation is 
the fact, that the worldline formulation of field theoretic correlators can be understood as 
the infinite string-tension limit of string-theory amplitudes, giving rise to a higher level of 
organization of the expressions. Important progress was made by generalizing the worldline 
techniques to effective-action computations [TJ [HJ - By now, the formalism has found a wide 
range of application, see, for instance, |H1 HOI HH [T21 [13 HH1 HZl CHI - The P ower of the 

worldline approach became particularly apparent in combination with numerical Monte- 
Carlo techniques [TTHl2()|: worldline numerics in the form of efficient and fast algorithms is 
now available for the computation of effective actions or Casimir energies for highly general 
forms of the external field. For instance, nonlocal phenomena and fluctuation-induced 
geometry-dependent interactions can now be computed reliably and straightforwardly |2*T) 

E2j. 

The majority of applications of the worldline approach remains perturbative in the cou- 
pling, even though formal all-order expressions can already be found in the early works j3] . 
The present work aims at exploiting these and related formal expressions for a direct eval- 
uation of all-order results by numerical Monte-Carlo means. We demonstrate the potential 
of worldline numerics for nonperturbative problems with the aid of a super-renormalizable 
model field theory involving a "charged" real scalar field <fi coupled to a "scalar photon" 
A. 2 Promoting the field formally to an iVf-component field, we study the system by 
means of an expansion for small Nf. Already at leading nontrivial order, this expansion 
corresponds to an infinite set of Feynman diagrams involving one open line or, alterna- 
tively, one closed <fi loop but arbitrarily many scalar photon fluctuations. We refer to the 
leading nontrivial order as the "quenched approximation". This approximation holds for 
arbitrarily large values of the coupling. Although it also includes infinitely many couplings 
to the external field, we expect that the approximation is most reliable for weak field am- 

1 Even though the method takes couplings to the external field into account to all-orders, the approxi- 
mations involved apply best to the case of weak fields, see below. 

2 The attributes in quotes serve only as an illustrative analogy, but should not be taken literally, since 
the model does not have continuous local nor global symmetries. Our model can also be viewed as a 
modified Wick-Cutkosky |28j model with imaginary coupling. It mimics as well the ubiquitous Yukawa 
coupling. 
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plitudes. As a first concrete application, we compute the nonperturbative effective action 
for the scalar photon, T[A], in the quenched approximation, concentrating on the induced 
interactions of soft photons. In this way, we obtain an effective action of Heisenberg-Euler 
type to leading order in Nf but to all orders in the coupling. 

The main purpose of the present work is to demonstrate that our method is capable 
of giving answers to nonperturbative problems with a computation from first principles. 
With regard to more realistic quantum field theories, our work serves as a feasibility study. 
Nevertheless, even in the simple model we observe new phenomena that may generalize 
to other theories as well. Of special relevance is the following result: whereas the pertur- 
bative expansion is ill-defined for a massless <fi field (similar, e.g., to massless QED), the 
nonperturbative effective action remains valid for massless fields, revealing the unexpected 
breakdown of perturbation theory (and not of the theory itself) in this limit. 

Our work is related to other nonperturbative investigations of field theoretic problems 
based on the worldline approach: in [21], scalar QED (ScQED) in quenched approximation 
was considered, and an instanton approximation of the path integral was used to compute 
all-order corrections to the Schwinger pair-production rate in a constant electric field. 
Furthermore, a comprehensive study of the worldline expressions for various quenched 
propagators has been performed in [23], employing a variational approach. Therein, the 
nonperturbatively interacting path integrals are approximated by optimizing a trial action 
according to a variational principle. In j2E], the worldline approach was used to deduce 
information about the nonperturbative quenched propagator beyond the well-controlled 
infrared (IR) limit. A combination of Monte-Carlo techniques with the worldline approach 
has also been used in j2Z] for the determination of bound-state properties from 4-point 
correlators; see also (2H1 for an approach to QCD bound states. Moreover, non perturbative 
worldline methods have been used to calculate the chiral anomalies pH].|rU)j. 

We start the presentation of our work in Sect.Elby rederiving the worldline expressions 
for the effective action and the nonperturbative propagator in quenched approximation. 
Even though the final results are known in the literature, we present their derivation here 
(with details in App. [A} in the modern language of generating functionals of correlation 
functions using functional integrals. In Sect. El we elucidate the numerical algorithm, 
concentrating on the new aspects introduced by the present work, and presenting results 
for the free theory as a useful example; another illustrative example is presented in App. 
C. In Sect. HJ we summarize our main numerical results for the self-interaction potential of 
the worldlines which lies at the heart of the nonperturbative dynamics. Section |S] employs 
this numerical data to derive our main result, the effective action of the scalar photon. We 
explain how the standard renormalization procedure applies to the present nonperturbative 
study. Section El summarizes our conclusions and gives an outlook. 



3 



2 Nonperturbative worldline methods in QFT 



Let us consider a specific Euclidean quantum field theoretic model in D dimensional space- 
time involving two interacting real scalar fields and A with bare Lagrangian 

£(0, A) = l -{d,<PY + \m 2 <? + \{d,Af - I h Atf. (1) 

This model meets our requirements in many respects: (i) it has a well-defined perturba- 
tive expansion, (ii) it will turn out to have also a well-defined nonperturbative expansion 
in the quenched approximation, (iii) the interaction ~ hAifi 2 is super-renormalizable in 
D = 4, since the coupling has positive mass dimension [h] = 1, and hence facilitates a 
simple analysis of the (nonperturbative) divergencies, (iv) the imaginary interaction imi- 
tates the interaction occurring in QED with A being the scalar analogue of the massless 
photon field; hence, A will be called the "scalar photon" in the following. We are well 
aware of the fact that this Lagrangian does not define a fully consistent QFT, since the 
potential is not bounded from below; however, contrary to the corresponding model with a 
real interaction (the Wick-Cutkosky model), the present model does not have propagator 
singularities which would signal the metastability of the model immediately and require 
a careful treatment. Nevertheless, the Euclidean effective action will turn out to exhibit 
nonzero imaginary parts which suggest a relation to decay widths of the metastable states. 
Note that in (21] it was argued that, despite the obvious non-hermiticity of the Hamiltonian 
of this model, the spectrum could be positive definite by virtue of a careful implementation 
of PT symmetry. 

Let us first concentrate on the effective action T[A] for soft scalar photons with momenta 
p 2 smaller than the scale of the massive field, p 2 <C m 2 . This effective action can be 
derived from the Schwinger functional, 

Z[J]=e w ^ = JvAV4>e-f* c ^ + f* jA 

= J VA det - 1 / 2 ( - d 2 + m 2 + ihA) e~ !* U^) 2 +L ™ 

VAe- s ^ (A)+ ^ jA , (2) 



where J = J d D x and S e R denotes an auxiliary effective action governing the quantum 
dynamics of the scalar photon (it should not be confused with the full quantum effective 
action T[A\ defined below). Here, we have integrated out the massive <fi field, resulting in 
a functional determinant that can be displayed in worldline form: 



S cS (A) = jl^A) 2 + l -\ndet{-d 2 + m 2 + ihA) 



-(dA) 2 - r ^ e~^T/ ihJ T drA( X (r))\ (3) 

2 v m / 2(Arr) D / 2 J T 1+D / 2 \ U [ ) 
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Here, the expression (. . . ) x denotes an expectation value with respect to an ensemble of 
worldlines x(t) with a Gaufiian velocity distribution, 



<•••>* 



x 



fVx< 



e~zfo dTr - 



(4) 



In a standard fashion, the effective action T[A] - the generating functional of all 1PI 
correlation functions of the scalar photons - is related to the Schwinger functional by a 
Legendre transformation, 

T[A] = -W[J] + / J A, (5) 



with J = J [A] being a functional of the "classical" field A (conjugate to the source J) 
defined by 

5W[J] 



A(x) 



SJ(x) 



= A[J\. 



(6) 



In this work, we will evaluate T[A] in the so-called quenched approximation. In a diagram- 
matic language, the quenched approximation for T[A] can be characterized by taking into 
account all diagrams with exactly one closed </> loop, but arbitrarily many virtual scalar 
photon fluctuations, see Fig. From a different viewpoint, the quenched approximation 
would correspond to the lowest nontrivial order in an N{ expansion if we promoted the 
field to have N{ components. 








• • • 



Figure 1: Diagrams included in the quenched approximation. 
More formally, let us introduce the field worldline current 



:= ih J dr5 {D \ 



x — x(t)), 



(7) 



such that the interaction of the scalar photon A with the <ft worldline can be expressed as 



J.J X J 



Then the quenched approximation consists in dropping all terms involving worldline cor- 
relators between two different j currents of the <p field, 0((ji . . .j2) Xl ,x 2 ) ~^ 0- Since each 
current couples to all orders to the classical field A, we expect this approximation to be 
good if the scalar photon field under consideration is weak compared to other mass scales 
of the system. Nevertheless, the quenched approximation does not restrict the value of the 
coupling h which therefore can be arbitrarily large. 
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In appendix [21 we derive the general worldline representation of the Schwinger func- 
tional and perform the Legendre transform yielding the effective action. In the quenched 
approximation, T[A] can be displayed in closed form (cf. Eq. (|75j) ). 

r nA L41 = / -id A) 2 - f T dT e - m2T (e^ f * iAj+ f* jA \ (9) 

QA[ l J x 2 { tl ) 2(4tt) d / 2 J T^/ 2 \ /*' { > 

where A = l/(— d 2 ) denotes the propagator of the scalar photon with coordinate space 
representation 

r(£=z) i 

A ( x i> x *) = -tm- \~.-~~\D-2 - ( 10 ) 



Fi — X2\ 



It is important to note that, in particular, all contributions from the scalar photon fluc- 
tuations are summarized in the exponential with a worldline current- current interaction 
~ f jAj. Re-expressed in terms of the worldline coordinates x(t) of the field fluctua- 
tions, this term acts like a self- interaction potential V[x] of the worldline, 



\D-2 



'0 Fl — X 2\ 

-vVf-r]. with. 7 :=^^r(^). (11) 



67T 

Here, we have introduced an effective coupling g for convenience. Therefore, a thorough 
analysis of the properties of this potential and its vacuum expectation value is mandatory 
for an understanding of the nonperturbative dynamics. This is one of the main aims of 
the present work. In addition to the effective action for the scalar photon, also proper 
vertices with external lines can be studied nonperturbatively on the worldline. In partic- 
ular, the computation of the propagator of the <fi field is straightforward and its worldline 
representation can be compactly written as 

Here, the expectation value (. . . )^ implies that we are dealing with open worldlines ranging 
from x\ to Xp. Again a proper understanding of the properties of the potential term is 
required. 

Furthermore, we would like to point out that the propagator can even be generalized 
to the case of arbitrarily many couplings to an external scalar photon background A by 
including a factor of e 19 ^ drA { x ( r )) j n the expectation value of Eq. ([12)1 . 

Let us conclude this section by translating the present results to scalar quantum elec- 
trodynamics (ScQED). In this case, the bare Lagrangian is given by 

AicQED = \DMW + m2 \<P\ 2 + \ F ^ F ^ (13) 
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where now denotes a complex scalar field and -D M [A] = <9 M + ieA^ is the covariant 
derivative. Most of our formulas can immediately be generalized to this case with the 
effective photon action given by 



l F F L_ 

4 ^ ^ (Aix) D l 2 



dT 



o 



jn+D/2 



C2 



2 fx ^M-^M^i^+Xc if^-M 



with the worldline current of the charged field 

j^x) = le 



r T 

/ rfri M (r)5 (1)) (x 



*(r)), 



(14) 



(15) 



and the photon propagator 

r(— ) 



1 + a 



Fi _ X 2 



D-2 ^ V 2 



a) 



(a?i - x 2 ) M (xi -ar 2 ) w 



D 



(16) 



with gauge parameter a. The interaction between the charged worldline current and the 
photon field then occurs in the form of a Wegner- Wilson loop, 



jef dx^A^x) 



(17) 



Here, the current-current interaction induced by the virtual photon exchanges inside the 
loop now implies a worldline self-interaction potential of the form 



with 

V{3cQED [x] 



dT\dr 2 



- I ■ A ■ 

2 J Jfi^livJv 



1 + a Xi ■ x 2 



•57T 



—r(^)v ScQEB [x] 



\Xl - x 2 



, • (4f- a) 



Xi ■ (x l - x 2 )(x l - x 2 ) ■ x 2 



\x 1 - x 2 



D 



(19) 

Note that the second term drops out completely in the Feynman gauge with a = 1. In fact, 
all a-dependent terms can be shown to correspond to a total derivative that automatically 
vanishes in the case of the photon effective action where all worldlines are closed loops jH] . 

The representation of the scalar propagator in Eq. (fT2*J) does also hold in the case of 
ScQED with V[x] replaced by VscQedN of Eq. (JTI| and #s c qed = ^byil^^). 

In the case of ordinary spinor QED, an additional complication arises from the cou- 
pling between fermionic spin and the electromagnetic field. Using a convenient spin-factor 
representation, the present formulation can be shown to hold with the simple insertion of 
the spin factor in all worldline expectation values, see |SSj- 



3 Worldline numerics 

In this section, we develop the discretization of the integral in the proper time. Since the 
discretization is done for an auxiliary variable, i.e., an integration parameter, it has the 
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advantage of preserving the continuum spacetime symmetries, i.e., the Lorentz (Euclidean) 
invariance. Also gauge and chiral symmetries would be preserved if present. This is 
the major difference between nonperturbative worldline numerics and other Monte-Carlo 
methods such as lattice field theory. In this sense, our method is much closer to analytic 
continuum field-theory techniques than to numerical lattice formulations. 



3.1 Algorithmic principles 

In this work, we evaluate the worldline integrals, corresponding to quantum mechanical 
path integrals, by means of Monte-Carlo algorithms, as proposed for perturbative ampli- 
tudes in PHI- For this, we approximate the infinite set of worldlines by a finite ensemble with 
til configurations. Furthermore, we discretize the propertime parameter of each worldline, 
such that the fluctuating worldline is specified by a set of N points per line (ppl), 

(...). = JP'--_fl'U J.£. ... with P[,] = e -w-, 
J VxP[x\ n L 

x{t) -> Xi e R D , i = 0,1,..., N. (20) 

In the case of the effective action, the worldlines form closed loops and we identify xq with 
Xjy; in the case of the propagator, where the initial and final point of the worldline are 
fixed, we use the definition x = x\ and x^+i = x-p- 

The finite number of worldlines introduces a statistical error that will be estimated 
automatically by the Monte-Carlo algorithm. The finite number of points per worldline N 
implies a systematic error that has to be analyzed straightforwardly by approaching the 
propertime continuum limit N — > oo. This is the most sumptuous point of the present 
investigation. 

An immediate realization of Eq. (|20|) requires to generate a different ensemble for each 
value of the propertime parameter T, occurring in the Gaufiian velocity distribution. This 
can elegantly be circumvented by introducing unit loops or lines y(t) 19 , 

y(t):=^=x(Tt), te[0,l], (21) 



such that the velocity distribution for y(t) becomes independent of T, 



drx 2 (r) = / dty\t). (22) 
o Jo 

Here, the dot always denotes a derivative with respect to the corresponding argument. As 
a consequence, one and the same ensemble of unit loops or lines can be used for all values 
of T, saving an enormous amount of CPU time. Of course, this change of variables also 
affects the remaining worldline integrand, for instance, the self-interaction potential, 

V[x] = [ dr x dr 2 = T ^ D/2 [ dt idt 



jo \xi -x 2 \ u Jo \yi -yi\ u 

T 3 - D / 2 v[y] (23) 
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where we have introduced the dimensionless self-interaction potential of the unit loops, 
v[y}. 3 For computations of the propagator, the variable transformation also affects the 
boundary points, inducing a further T dependence, xp = VTi/f, x\ = y/Tyi, that can 
nevertheless be handled straightforwardly. 

In the context of the effective-action calculation involving closed worldlines, one final 
convenient step for the algorithm consists in using the center-of-mass decomposition of the 
worldlines, x(Tt) = x CM + VTy(t), such that the unit loops satisfy f Q dty(t) = 0. This 
implies that the ensemble average becomes 

(F[x]) x = J d D x CM (F{xcm + Vfy]) y . (24) 

This allows to study the effective action density, i.e., the effective Lagrangian, locally, 
with the center-of-mass integration corresponding to the usual relation between action and 
Lagrangian, r = J d D xcM£- 

For the Monte-Carlo generation of the worldline ensembles, there are various efficient 
techniques available that do not require to spend computer time on dummy thermalization 
sweeps, as it would be the case for standard heat-bath algorithms. Fast algorithms are, for 
instance, the Fourier loop algorithm, "f loops" , or explicit diagonalization of the velocity 
distribution, "v loops" , see [22] • For most computations in the present work, we use a new 
algorithm that is based on a point doubling procedure ("d loops" and "d lines") that is 
described in appendix [Bj 

Beyond these general algorithmic aspects, we encounter a particular numerical problem 
in the present case upon discretizing the self-interaction potential ()23|) : the integrand for 
coinciding points y± = y% is ill-defined. The possibly arising divergence can be viewed 
as a "self-energy" of the worldline. We regularize the discretized sum by taking out the 
coincident points by hand, 

"M-^E l *p-2 » V - (25) 

A singular behavior of the original integral then manifests itself in a strong N dependence 
of v[y] which we will carefully control. Since the distance of the nearest-neighbor points 
decreases for increasing N, \yi — ~ 1/y/N, such strong TV dependence signals a short- 
distance UV singularity that will be subject to renormalization of physical parameters. 
Alternatively, we could include the coincident points in the above sum by shifting the 
denominator, \yi~ yj\ — > \yi — yj \ +e, and study the limiting behavior with e — ► 0. However, 
since the latter limit will always interfere numerically with the large- iV limit (which has to 
be taken anyway), the first proposed method can be handled more straightforwardly. 

3 Incidentally, we note that the same transformation implies for the potential in ScQED that VscQED [x] = 
T 2 ~ D / 2 vs c QEn[y}', therefore, Vs c qed is already dimensionless in D = 4 and the unit-loop transformed 
exponential of the potential becomes independent of T . 
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3.2 Free theory 

As a first useful application, let us study the worldline representation for the free theory in 
the limit of vanishing interaction, e.g., h — > for the purely scalar system. An instructive 
example is given by the kinetic action for the worldlines, S w , characterizing the Gaufiian 
probability distribution of the worldline ensemble, 

P[x] = ATexp(-S w [x]), (26) 

where M is the normalization, ensuring J T>x P[x]—1. Discretizing the free worldline action 
S w (introduced in Eq. (jl))), as outlined above for the case of the propagator, we obtain 

1 f T N 4- 1 N+1 

S, = - drx 2 (r) -> — Z- ~ x^)\ (27) 

^° i=i 

with and xn+i = xp the corresponding initial and final points of the propagation. 

For the propagator, we would finally have to integrate over the propertime T. However, 
we will be satisfied here with the propertime integrand and set T to some fixed value in 
the following. 

The ensemble average and the value of the root mean square (RMS) of the free worldline 
action can be calculated analytically. For this, let us consider the auxiliary integral 

dx x ■ ■ -dx N exp{ — 22 1 Ar 2~ — }' ( 28 ) 

i=i 

where At = T/(N + 1), and we have suppressed the Lorentz indices of the a^'s which 
are points in D dimensional spacetime. The parameter k is arbitrary here which will be 
exploited below. 

The integral is Gaufiian and can be done exactly. We rewrite the integral as 

/kAt 1 
dx\ ■ ■ ■ dxivexp{ — —-^{x\-\-x'^ — 2xiXi—2x^x^-\-XiAijXj)}, (29) 

where A is a symmetric matrix, 

A {j = 2Sij - S i+1J - 6 i>j+1 , i, j = 1, . . . , N. (30) 
Introducing the auxiliary current rji, 

T)i = -2xi5n - 2x F 5 iN , (31) 

the integral reads 



AttAt jvd 1 
I N (x F ,xi,T,K) = { ) 2 n exp 

K ^/detp) 



K f]i(A 1 )ijf]j K ^2 _^ x 2 ^ 



AAt 4 4Ar 



(32) 
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N 


2T 


Xp 


( 'S'w ) numerical 


((S'w) theory 


AS W /S W (%) 


31 


0.5 


(10,0,10,0) 


161.895 


162 


-0.06 


1023 


0.5 


(0,0,0,0) 


2046.38 


2046 


0.02 


1023 


0.5 


(10,0,0,0) 


2096.00 


2096 


0.00 


4095 


0.5 


(0,0,0,0) 


8190.15 


8190 


0.001 


4095 


0.5 


(10,0,0,0) 


8238.70 


8240 


-0.02 


4095 


0.5 


(100,0,0,0) 


13189.4 


13190 


-0.005 



Table 1: Average theoretical and numerical actions. All numerical calculations were done 
with n L = 10 4 worldlines. We have set D = 4 and x\ = (0, 0, 0, 0). 

with 

Vi(^~ 1 )ijVj = 4(xi(A- 1 ) 11 + 2x 1 x F (A~ l ) 1N + A^ N x 2 F ). (33) 
The ensemble average of the free action can now be calculated from the auxiliary integral, 



(S w ) = M [ dx\ ■ ■ ■ dx N S w exp{— S w \ — N [ dxi ■ ■ ■ dx N — — exp{— kS w } 
J J dn 

The integral actually agrees with 1^ of Eq. (f3~2"J) above, therefore 

ND 

.-, —d ( 47rAr\ 2 1 r k 



■ (34) 

K = l 



K=l 



= \nD + ^(x f - Xi ) 2 . (35) 

Notice that the result corresponds to the classical action plus one half times the number 
of degrees of freedom. 

In tabled we show a comparison between the exact results and those of the numerical 
calculation using our Monte-Carlo code for different parameters. The numerical result is 
always very satisfactory with a precision better than 0.1 % in all cases. This illustrates 
the power of the Monte-Carlo method. We underline a fact that will be of the outmost 
importance in the more involved non perturbative calculation below, namely that the 
Monte-Carlo method also provides information about the probability distribution of the 
free action itself. For instance, the probability of finding a configuration with worldline 
action S' w is given by 

V(S' W ) = AT J dx x ■ ■ ■ dx N 5(S W - S'J exp{-S w }, (36) 

which can be evaluated analytically, yielding 



ND 



V(S'J = Af(S' w - S dass )- 1+ ^e~^~ s ^ 6(S' W - S class ), (37) 



11 



where Sdass = (l/4T)(xp — x\) 2 is the classical action. In Fig. El we compare the action 
distribution calculated with our Monte-Carlo code to the analytical result. The agreement 
with the theoretical calculation is excellent. 4 Using the probability distribution V(S' W ), it 




Figure 2: Action distribution from the Monte-Carlo method (histogram) in comparison 
with the analytical calculation (continuous line) for N = 31 and x-p = (10, 0, 10, 0). 

is straightforward to calculate the ensemble average, the RMS, and higher moments of the 
free action. The ensemble average immediately verifies Eq. (jSSJ), whereas the RMS results 
in 

* 2 (S W ) = (Si) - (S w ) 2 = (38) 

As N increases, the ensemble average S w (xf,xi) grows as well, and the distribution of S w 
broadens. For large N, we observe from 

- ~A=r ( 39 ) 
(S w ) yfND/2 

that the relative width of the distribution becomes narrower. 



4 Self-interaction potential in the scalar model 

In this section, we perform a detailed numerical analysis of the self-interaction potential, 
occurring in Eq. (jlljh In particular, we consider its dimensionless form v[y] defined in 
Eq. (j2HJ) in its discretized regularized version, 

^Sra- (40) 

4 We have chosen N small to exhibit more clearly the non Gaufiian character of the distribution 
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as discussed in Eq. ()25p . In the following, we concentrate on D = 4 dimensional spacetime. 
Other dimensions can be treated in the same way. 

The average value of the self- interaction potential v[y] as a function of the number 
of points per loop can be seen in Fig. Altogether with a logarithmic fit. Obviously, the 
self-interaction potential diverges logarithmically with N, 

(v) y = a + b\n(N), (41) 

with a ~ 0.363 and b ~ 0.341. By contrast, the RMS is finite, as is visible in Fig. 0] We 
have furthermore collected numerical evidence that all the higher-order moments are also 
finite. 




6 8 10 12 14 



log 2 (N) 



Figure 3: Average value of the self- interaction potential as a function of the logarithmic of 
the number of points per loop (iV). 

In order to evaluate ensemble averages, it will become useful to study the distribution 
of the self-interaction potential defined by 

P( v ) =ft Jvy5(v- v[y]) exp(-S w [y] -v[y}), (42) 

where the normalization is fixed such that J dvP(v) = 1. We depict the numerical result 
for this distribution in Fig.Elwith binned values of the potential for two worldline ensembles 
with different iV ppl. The distributions are shifted along the x axis by the logarithmically 
divergent ensemble average of v , but are not modified otherwise. As a result, all the distri- 
butions are equal within the numerical accuracy for different values of N, demonstrating 
that all other properties of the distribution are devoid of further divergencies. 

This remarkable property reflects the fact that we are dealing with a super-renormalizable 
interaction. In a diagrammatic language, only a finite subset of diagrams is superficially 
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Figure 4: RMS of the scalar to scalar potential as a function of the number of points per 
loop (N). 



divergent. In the present scalar model, this finite subset consists of the one-loop contri- 
bution to the 2-point function (mass operator), the one- and two- loop scalar photon 
tadpole and the photon 2-point function ("vacuum polarization"). Once these diagrams 
are controlled by adjusting the counterterms and renormalizing the corresponding physical 
parameters, all higher-loop diagrams do not contribute independently to further renormal- 
ization; divergencies in subdiagrams will be canceled by the already adjusted counterterms. 
Evidently, our nonperturbative evaluation of the infinite set of diagrams contributing to 
the quenched approximation will also benefit from the simplified renormalization in this 
particular scalar theory. 

Let us show that all the characteristics of the average potential can be understood in 
terms of the distribution of chord distances and that this distribution gives an intuitive 



Figure 5: Distribution of potentials for N = 2 9 (continuous line) and iV = 2 12 (dashed 
line). The average values of the distributions have been shifted to lie on top of each other. 
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picture of the problem. For a large number of points per loop we can introduce the chord- 
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Figure 6: Chord-distance distribution between pair of points in paths for several values of 
points per loop, as marked. 

distance distribution, 




where = |^ — yj\ abbreviates the chord distances and p(l) is the probability of having 
two points separated by a distance I. 

In Fig. IHl we show the chord- distance distribution in a double-log plot for N = 2 5 ,2 12 , 
and 2 14 . From the figure, we can see that very large chord distances are exponentially 
suppressed. For very small values of the distance, / <C l/\fN, much less than a typical 
random step, we see that the probability goes like which is simply a measure of the 
available volume in D dimensions. This behavior is a direct consequence of the world- 
line discretization, and the corresponding branch vanishes in the continuum limit. For 
intermediate distances the distribution goes like Z 1 which implies that the path spans a 
2 dimensional surface, characteristic of a random walk. The contribution to the average 
potential of the very low chord- distance region is finite, in 4 dimensions. However the 
intermediate region gives a contribution of order 

f dlp(l)^~-\n(l/VN), (44) 

which diverges in the continuum limit and gives rise to divergencies in the average potential 
and, consequently, in the effective action. The divergent part of the average potential 
carries information about the topological properties of the paths which contribute to the 
average. 
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5 Effective action for the scalar photon 



Let us now exploit our knowledge about the worldline self-interaction potential in the scalar 
model in order to compute the nonperturbative effective action T[A] of the scalar photon in 
quenched approximation in D = 4 spacetime dimensions. In this section, we demonstrate 
how the renormalization procedure can be carried out to all orders, using this example. 
For this, we need to consider the properties of the potential V[x] for closed worldlines, or 
its dimensionless counterpart for unit loops, v[y], cf. Eq. 

In addition to the quenched approximation, we confine ourselves to slowly varying pho- 
ton fields with frequencies smaller than the mass scale, uj 2 <C m 2 , for which A can be 
assumed constant in spacetime. Compared to QED, these assumptions are reminiscent to 
those leading to the famous Heisenberg-Euler effective action. Notice that this approxima- 
tion is consistent with the quenched approximation. 

Employing the representation (JHJ) together with Eq. (fTT|) and the unit-loop transforma- 
tion of Eqs. ()21j) - ()23j) . the effective action in these limits reads in D = 4, 



where we have taken the photon field dependence out of the worldline average, owing to 
the A = const, assumption. 

We make use of the potential distribution, as introduced in the previous section, since 
it is sufficient for evaluating the nonperturbative effective action of the scalar photon and 
performing the renormalization program. Let us illustrate this statement with the aid of a 
simplified fit to the potential distribution; a better but more complex fit will be presented 
below. As a simplified parameterization of the potential distribution, let us consider 

P ^ = ^f-Tu ( y ~ - wo)) 9(v - v ), (46) 

1 [a + 1) 

where a, f3 and v denote fit parameters and 9 is the step function. This distribution is 
already normalized, J dvP(v) = 1. Numerically, we find that a and (3 are finite, a ~ 0.79 
and (3 ~ 13.2. Only vq diverges with the number of points per loop, which can be related 
to the divergence of the potential ensemble average, since 

< v > y = J dvv P(v) = v + ^p- (47) 

With regard to our result of Eq. f)41j) . this implies that 

1 + a 

v = blnN + a — ~ 0.34 In iV + 0.23 . (48) 

P 

Most importantly, the distribution also provides us with a numerical result for the expo- 
nential of the potential, carrying the nonperturbative information about the infinite set of 
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quenched diagrams, 

( e -flT«M) i/ = J dvP(v)e- gTv 



) e-^ = F M) (^T)e-^ . (49) 



^)(^)=(^-l) • (50) 



For later convenience, we have introduced the auxiliary function 

Inserting this result into Eq. (|45j) . we obtain the numerical estimate for the unrenormalized 
effective action, 

r Q A[A] =~J d*xJ*^e-^^F {aP) {gT)e-^, g = ^. (51) 

As for the scalar-photon operators, we demand the renormalization conditions (indepen- 
dently of the approximations) 



5T[A = 0] 5 2 T[A = 0] 



<L4(x) ' <L4(-p)<L4(p) 



0. (52) 



The first condition implements the "no-tadpole" renormalization prescription, whereas the 
second condition defines the scalar photon to be massless. In general, a third condition 
fixing the wave function renormalization of the scalar photon Za has to be implemented; 
however, it is a peculiarity of the present model that the renormalization shift of this 
Za factor remains zero 5Za = 0. 5 Incidentally, the same is true for the wave function 
renormalization of the field. 

These renormalization conditions fix the corresponding counterterms, resulting in the 
following (partially renormalized) effective action: 

Tqa[A] =~J d*x £ ^ e- 2T (e^ - 1 - ihAT + <*^2!) F^gT) e~^, 

(53) 

where we have also subtracted the zero-point energy, enforcing T[A = 0] =0. Note that 
these renormalization conditions have rendered the propertime integral finite. This is, in 
particular, independent of the loop order in the quenched approximation, since the coun- 
terterms do not receive contributions from higher-loop orders due to super-renormalizability. 6 . 

The remaining divergence contained in vq can finally be removed by mass renormal- 
ization. In principle, this is done by imposing a renormalization condition for the 2-point 

5 This is reminiscent of the same relation of 4 theory in D = 4 to one-loop order; because of the 
super-renormalizability of the present model, this relation holds to all orders here. 
6 This point is expected to be different in renormalizable theories, such as QED. 
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function of the field at zero momentum. This is straightforwardly possible within our 
approach but requires an explicit computation of the 2-point function. However, we adopt 
a simpler prescription here for which it suffices to consider merely the scalar photon action: 
we demand that the propertime integrand should fall off exponentially for large T with a 
width given by the renormalized mass, 



^ e~ m2T (j hAT - 1 - ihAT + ^2!) F (aP) {gT) e~^° 



~/(T)e- m ^, (54) 



where the prefactor f(T) has a weaker T dependence than an exponential function. This 
renormalization condition also fixes the finite part of the mass renormalization according 

to 

2 h 2 2 ( 1 + cA 
m R = m +——v = m +g [bhiN + a — . (55) 

Q1T Z \ fj J 

We end up with the fully renormalized effective action for soft scalar photons in the 
quenched approximation, 

W4 = / jf g ^ - 1 - **r + (£r 

, (56) 

where the form of the function Ft a p\ (x) is given in Eq. (JSTij) and holds for our simplified fit 
(JUJ). 

A more accurate parameterization of the potential distribution can be obtained by the 
fit function ^ 

P(v) = pi exp{-— -J (log(v - v ) - p 3 ) 2 } 0(v - v ), (57) 

a plot of which is shown in Fig. [Altogether with the potential distribution for N = 2 14 . The 
Pi as well as Vq are the fit parameters. Again, we observe that only v carries a logarithmic 
divergence, whereas the Pi are finite numbers. From the fit we obtain p 2 — 0.24, p 3 ~ 0.60, 
and pi ~ 0.89 fixes the overall normalization, vq is related to the average potential 

Vo =< v > — L e P3+3/2 P l _ 0.34 ln(JV) - 0.76 . (58) 
v 71 

This improved fit results in a better numerical estimate of the potential average, 



o 



= : e~° Tv ° F^gT), (59) 

which is of the same form as for the simplified fit in Eq. (|4*9~|) . but involves a more complex 
function F( p .)(x) defined by a parameter integral. The mass renormalization proceeds as 
above, = m 2 + h 2 / (87r 2 ) v , since F^(gT) decays weaker than exponentially for large T, 



18 



3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 £ 

v(y) 



Figure 7: Probability distribution of the self-interaction potential (histogram) and a fit to 
a modified log Gaufiian (continuous line). 

see App. |D| Therefore, replacing F( a p)(x) in Eq. (j^Bj) by F( Pi )(x), we obtain our improved 
result for the effective action for the scalar photon. 

At this point, it is important to stress that both fits do not only differ with respect 
to their quality of describing the data, but also with respect to their analytic properties. 
This becomes particularly apparent from the corresponding perturbative expansions. For 
instance, the simple fit leads to F^ a ^(gT), Eq. ()50|) . which is completely analytic in g and 
hence can be expanded in a well-converging Taylor series, as displayed in Eq. (jl(J4|) in 
App. |Dl By contrast, the resulting auxiliary function F^(gT) from the better fit gives 
rise to an asymptotic series upon expansion for g — > + , as shown in Eq. (jlOHj) . 

Within the language of diagrammatic perturbation theory, the implications from the 
two fits are different. To see this, let us note that the ratio of successive loop contributions 
to a scalar photon proper vertex (obtained by expanding Tqa[^4] in terms of A) is pro- 
portional to the ratio of the corresponding expansion coefficients of the auxiliary function, 



r(") I 

1 QA,RU+1 loop 



r(") I 

1 QA,RU loop 



-{l + n) 



l.g. in I 



F®(gT = 0) 
F (i-i)( g T = 0) 



QA,R 



QA,R 



[A = 0] 



5 A... 5 A 



(60) 



where F^ l \gT = 0) denotes the Zth expansion coefficient of the corresponding auxiliary 
function (see, Eqs. (jl03J) and ()104J0 . and we have concentrated on the leading growth (l.g.) 
of the ratio for large /. The prefactor of (/ + n) arises from the propertime integration and 
can be traced back to the fact that the self-interaction potential V is dimensionful in the 
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scalar model. 7 Now, for the two different fits, this coefficient ratio reads 



r Q \ R | m ioo P f if for F( a p) 

r (n) , ^ X 1 e Pia+l/2)- P 3( 1 + n ) f p ■ l Di J 

1 QA.RU loop l g inl < V i J yp l ) 

This implies that the proliferation of diagrams with increasing loop order translates to an 
increase of the expansion coefficient which is very different for the two fits. The simple fit, 
Eq. (|46p. leads to a mild increase of the expansion coefficients. Though the resulting per- 
turbative series would not be convergent but asymptotic, the series would be immediately 
Borel summable. 8 For instance, for the case of ScQED, it has recently been conjectured 
that the perturbation series even converges for the on-shell renormalized QED n-photon 
amplitudes in the quenched approximation [34J. 

For the better fit, Eq. ([57)1 . the coefficient ratio increases strongly, implying that a 
resummation requires generalized resummation techniques. Since the leading-growth be- 
havior of the coefficients from the loop expansion are influenced by the combinatorics of 
diagrams, the conclusion for this perturbative combinatorics would be very different for 
the two fits. The lesson to be learned is that drawing a reliable conclusion about the 
analyticity structure of the perturbative expansion from the numerical estimate is very 
difficult. Different fits of similar quality can lead to very different analytic structures. 
However, since we have the fully resummed finite result numerically at our disposal, an 
analysis of the analytic structure may be viewed as a rather academic question. Notice 
also that the large gT expansion is controlled by the small v behavior, i.e., it is related to 
the approach to the "classical" region where v ~ f m m- This classical region corresponds to 
a subset of worldlines that are close to the minima of the action and thus resemble classical 
trajectories. 

Let us finally point out one remarkable conclusion that can be drawn form the effective 
action (|5fij) . To one-loop order, the small- mass or large-photon- field limit of the effective 
action is of the form 

i r h a 

TuiooM] I ha >M ^ ~-— 2 / d 4 x (hA) 2 In — , (62) 

exhibiting a logarithmic increase, as is familiar from the Heisenberg-Euler QED effective 
action or the Coleman- Weinberg potential. This form implies that the massless <fr field limit 
is not well defined for our scalar model; technically speaking, the T integral is divergent at 
the upper bound in the massless limit. This becomes different, once the infinite number of 
diagrams of the quenched approximation is taken into account. These induce the F^ ^(gT) 
function which vanishes for larger T, F^ ^(gT — > oo) — > 0. This property renders the T 
integrand finite even in the massless <fi field limit. For instance, for the simple fit, the 

7 This is different in ScQED where the potential is dimensionless in D — 4; hence, this factor ~ I would 
be absent. 

8 In the somewhat artificial large- n (large photon number) limit, i.e., u>l, the series is indeed absolutely 
convergent for small coupling, being of the type of a geometric series. 
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massless limit of Tqa[A] can be evaluated analytically. Here, we just cite the result for the 
limit A/h <C 1, for which the quenched approximation is expected to be more reliable, 

r^A]L^-- l ~ r ^l~^J a J d>x ( hA ~> 2 g)*[i + o(Wk))]. (63) 

The overall sign is, of course, still negative since a ~ 0.79. We conclude that a completely 
massless scalar model of the present type can be formulated consistently in the quenched 
approximation; it is only the perturbative expansion of this model that breaks down in the 
massless limit. 

6 Conclusions and Outlook 

The feasibility of a nonperturbative analysis in quantum field theory based on the worldline 
formalism has been demonstrated here. The approach combines the efficiency of worldline 
representations, i.e., having infinitely many Feynman diagrams in one worldline expression, 
with powerful numerical algorithms for the worldline integrals. As a major advantage 
of this approach, all continuous spacetime symmetries remain unaffected by numerical 
discretization procedures. 

In the present exploratory study, we have considered a scalar model with cubic inter- 
action in four dimensions in the quenched approximation which is the leading non-trivial 
order in a small- Nf expansion. This approximation is expected to hold for small external 
field amplitudes but arbitrary values of the coupling. The model mimics scalar QED or 
Yukawa theory and is super-renormalizable though not finite. Therefore the divergence 
structure is rich enough to allow for a detailed study of the main difficulty of many field 
theoretic methods: understanding the origin and taming of divergencies. The combination 
of the propertime representation and the worldline discretization provides for a sufficient 
regularization, facilitating the identification of divergencies. Subsequently, the divergen- 
cies can be removed by fixing physical parameters by virtue of renormalization conditions 
on a nonperturbative level. It is particularly such a nonperturbative analysis of quantum 
field theories which sheds light on the fundamental problem whether a perturbative for- 
mulation gives an exactly correct account of divergencies jHHj- For practical purposes, our 
divergence analysis may be helpful for technically similar problems, e.g., arising in many 
effective theories for polymer physics 39J. 

For the case of the effective action for soft scalar photons, we have performed this pro- 
gram explicitly. The result is a finite propertime parameter integral of Heisenberg-Euler 
type, containing all-order corrections in the coupling. These all-order contributions par- 
ticularly modify the large-propertime behavior, in turn affecting the small-mass limit. We 
observe that the zero-mass limit surprisingly exists in the nonperturbative result, whereas 
perturbation theory is ill-defined in this limit. Hence, we conclude that the large loga- 
rithms for small masses from all quenched diagrams can be summed up to a finite result. 
Whether this feature persists beyond the quenched approximation remains an open issue. 
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At this point, a more detailed discussion of the interplay between the various parameters 
and approximation methods is in order. In the scalar model, we have three massive scales, 
implying two dimensionless parameters, e.g., hA/m 2 and hA/h 2 = A/h. In the limit 
A/h ^> 1, radiative scalar photon exchanges are suppressed and the perturbative loop 
expansion is expected to become reliable, almost independently of the value of hA/m 2 ; of 
course, for hA/m 2 3> 1, infinitely many couplings to the external field have to be taken 
into account. For A/h <C 1, the systems is in the nonperturbative regime. In the additional 
limit of hA/m 2 <C 1, higher <fi loops are suppressed by the mass threshold in comparison 
to radiative photon exchanges; hence, the quenched approximation is reliable in this limit. 
Whether this suppression of <fi loops persists also in the massless limit of the scalar model 
is an open question. The fact that the massless limit exists in the quenched approximation 
may be taken as an indication for an affirmative answer. At least, there is no obvious 
reason for additional virtual <fi loops to exhibit a problematic infrared behavior; on the 
contrary, virtual momenta flowing into the virtual loops rather serve as additional infrared 
regulators. In this sense, the quenched approximation with one virtual <ft 1°°P an d soft 
zero-momentum external photons represents the case of highest sensitivity to the massless 
limit. To summarize, we expect the finiteness of the massless limit in the scalar model to 
persist beyond the quenched approximation; whether further <ft loops give quantitatively 
sizable contributions to this result has to be checked by explicit computation. For this, the 
necessary inclusion of additional internal matter loops into our formalism appears to be 
feasible with our present technology. We would finally like to stress that the massless limit 
of the scalar model may be a peculiarity induced by the existence of a massive coupling, 
and it might not translate to a similar property of renormalizable models such as QED 
or Yukawa interactions. For instance in spinor QED, an interacting massless limit of the 
soft-photon action is inhibited either by triviality or by chiral symmetry breaking [35J. 

From a technical perspective, we have developed further the algorithmic tools of world- 
line numerics. Beyond those techniques that so far have successfully been applied to 
perturbative computations, the use of probability distributions renders the algorithms 
extremely efficient. For instance, the effective action, i.e., the generating functional of 
infinitely many soft photon amplitudes, has been derived from one single distribution (of 
the self- interaction potential). A further illustrative example for the harmonic oscillator 
is given in App. O We are convinced that the method will find a variety of applications, 
given the large number of problems that are being tackled by propertime methods. 

Our work paves the way for many generalizations: for instance, both the worldline 
formalism as well as our worldline numerics are well suited for a study of propagators in 
coordinate space, where the scale set by the propagation distance represents one further 
technical challenge. The study of propagators appears mandatory for an analysis of renor- 
malizable theories such as (scalar) QED in D = 4. In these cases, our simplified mass 
renormalization will no longer be sufficient and the renormalized mass needs to be fixed 
with the aid of the propagator's correlation length. 

Moreover, a generalization to gauge interactions will require a more careful treatment 
as far as the discretization is concerned. Our present method if applied to a gauge theory 
with a charged field would not satisfy current conservation which is essential for main- 
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taining gauge invariance. More refined algorithms are required and are a subject of current 
development. 
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A Derivation of the effective action 

Here, we derive the worldline representation of the Schwinger functional and the effective 
action for the scalar photon field T[A], Eq. (J2J). We start by citing the Schwinger functional 
as given in Eq. (J2J) with the insertion of the worldline form of S e s of Eq. (JHJ), 



Z[J] 



,W[J] 



VA exp 



2(<W + 



+ / J A 



(64) 



where we have used the worldline current representation of Eq. (JZJ), and introduced the 
abbreviation 



1 



T ' 2(4tt) d / 2 J T x + d / 



dT 



-m 2 T 



(65) 



Let us first expand the worldline exponential, such that the A integral becomes Gaufiian 
at each order of the Taylor sum, 



Z[J) 



tv. 



Ti 



£ 



00 

Y- 

1,1 I 



VA e~ I* ^ 2 ef* Mh+-+3n+J) 



X\...X r , 



C2 



1 f K J AK J 

2 Ji I ^n'- ll ^n 



X\.. .X r , 



(66) 



where we have introduced the generating functional of connected Green's functions W[J] 
in the last line, and have defined the sum of currents 



(67) 



23 



The kernel A denotes again the scalar photon propagator as given in Eq. (|1U|). Note that 
Eq. ()66|) is an exact representation of the Schwinger functional without any approximation 
so far. The classical field A (conjugate to the source J) is obtained from W[J] by functional 
differentiation (cf. Eq. (jHJ)), 

5W[J] 

A ~ ~JT 

= AJ + e-^'^i/ /Aj2j ie y* K Z AK i\ . (68) 

n U ' ^ T i- T « \ i=l / xl ... Xn 

A partial inversion of this equation gives, 

J = A- 1 A-e- w J2^ I (X>e^^ A ^) , (69) 

with A -1 = — d 2 . Note that this is an implicit definition of J = J [A], since K% in the 
exponent on the right-hand side still depends on J. Nevertheless, these definitions can be 
used to give an implicit definition of the exact effective action for the scalar photons T[A] 
(cf. Eq. ©), 

T[A] = -W[J[A}} + [ AJ[A]. (70) 

J X 

Explicit representations can be found in the quenched approximation that consists of drop- 
ping all terms involving worldline correlators of different worldline currents, (ji . . . jj) Xi ... Xj — » 
0, i 7^ j. In this approximation, the Schwinger functional Eq. (jfifij) can be written as 



w QA [j] = e ±J x JAjJ2^ J (e\LhNi+S x JN^ J ^ \i,.j„±,„-i,j± J .^ 



exp 



(71) 



The corresponding representation of the classical field in quenched approximation boils 
down to 

j) • (72) 



A = AJ + j /a, A 



Note that, even in the quenched approximation, this equation cannot be inverted explicitly 
for the current, 



J = A~ 1 A- / / e U^+/«JAj A . (73) 



T 



We take this equation as the implicit definition of J = J[A]. By virtue of Eq. (fTUj). this 
allows to express the effective action Tqa[^4] intermediately in terms of J[A], 

T QA [A} = \J JAJ + J^^ A ^ JA ^(JAj - 1))^. (74) 
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Now, J [A] of Eq. (J75j) can be iteratively inserted into Eq. (J73J). As a consequence, the 
iteration terminates after the first order by means of the quenched approximation, since 
the higher orders involve multi-worldline current correlators. The result can be displayed 
in closed form, 

I*M = I \(^f - ^ [ ^( e iU^U*), (75, 
which serves as the starting point of our investigation in Eq. (JHJ). 



B Worldline generation: "d loops" 

Here we describe a new algorithm for generating open or closed worldlines ab initio, which 
is used for most numerical studies in the present work. The algorithm generates worldlines 
by a suitable doubling of points, starting from a small number (e.g., from 1 point); hence 
we call the resulting worldlines "d loops" (closed) or "d lines" (open). The algorithm is 
reminiscent of a heat-bath algorithm, as was used in earlier worldline numerical studies 
[TH] . but does not require dummy thermalization sweeps. 

Let us first consider closed "d loops". By appropriate rescaling, the Gaufiian velocity 
distribution can always be brought to unit loop form, 

P[ y ] = e -Uo^(^) 2 . (76) 

Discretizing the t derivative by an asymmetric nearest-neighbor form, dy/dt = (yi — 
2/i_i)/At, with i — 1, . . . , N , the distribution becomes 

PnM -> e-^ESifo-w-O 2 . (77) 
The probability for the location of the z'th point is 

py = e -^[2^ 2 -2j /l (j /l+1 +y l _ 1 )]... _ e -^2[j/ l -i(y l+1 +j /l _ 1 )] 2 ...^ 

such that i/i is in a GauBian heat bath of the algebraic mean of its nearest neighbors. The 
crucial point now is that we can add another N points in between with the following 
probability distribution, 

Vyiyi — e -^[2(y l ,-\{y l +y l - 1 )) 2 +2(y ll+1 -\{y l+1 +y l )) 2 }...^ ^ 

such that the new total probability distribution for yi becomes 

Pvi - = e- 2J ^ 2 ^^' + ^ 2 -. (80) 

Note that the probability distribution for the doubled points yy in Eq. (f?H|) has been 
adjusted such that the probability distribution for the point yi has lost its dependence on 
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y i+ i and yi-i, Hi only feels the heat bath of the new points yy and yy+i. Correspondingly, 
the full distribution for the iV + iV = 2N = Ni points reads: 

PnAv] -> e-^-^^-^- l] \ (81) 

where we have redefined yi — > y 2 i, yy — > y2i-i- 

Now the recipe for the "d loop" algorithm generating closed worldlines with N ppl is 
obvious: 

(1) Begin with one arbitrary point N = 1, y^. 

(2) Create an N\ = 2 loop, i.e., add a point yjv/2 that is distributed in the heat bath of 
yN with 

e ~^r 2 (yN/2-yN) 2 _ ^g2) 

(3) Iterate this procedure, creating an iV fc = 2 fc ppl loop by adding 2 fc_1 points y q N/N k , 
q = 1,3, . . . , N k — 1 with distribution 

(4) Terminate the procedure if N)~ has reached = N for unit loops with iV ppl. 

(5) For an ensemble with common center of mass, shift each whole loop accordingly. 

Open worldlines ("d lines") can be generated analogously by starting with the end 
points xi = yo and xf = yjv+i and filling N new points in between by the same doubling 
procedure. 

In comparison with existing worldline generator algorithms such as the "f loop" and 
"v loop" algorithms of [22] , all algorithms share the property that each generated random 
number is used for a worldline; in particular, no dummy thermalization sweeps are neces- 
sary. Similarly to the "f loop" algorithm, "d loops" produce worldlines with N = 2 l ppl, 
/ G N, whereas the "v loop" algorithm can produce any N G N. "f loops" by Fourier trans- 
formation discretize the derivative x implicitly such that the operator d 2 /dr 2 is diagonal 
in Fourier space, whereas "v loops" and "d loops" use the asymmetric nearest-neighbor 
discretization x ~ Xj — x^\. The later has proved to approach the continuum limit in prob- 
lems with scalar interactions faster than the Fourier decomposition [22] • in comparison to 
"v loops" , the present "d loop" algorithm requires less algebraic operations and hence is 
faster by some ten percent. 

C Harmonic oscillator 

As an instructive example of the use of the Monte-Carlo method involving distribution 
fits, let us study the harmonic oscillator. The propagator for a harmonic oscillator can be 
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written in the form 



G{T,R) 



x(T)=x f 



x(0)=xj 



Vxexp{- 



dr x + 00 x }, R = \xf — x\\. 



34) 



In the language of quantum field theory, this corresponds to the heat kernel of a D — 1 
dimensional scalar field theory, coupling to a constant "potential" uo 2 . In the notation of 
worldline ensemble averages, this propagator can be written as 9 



G(T, R) = ( exp 



7 2 2 

aruj x 



(85) 



Introducing the T scaling as outlined in Eq.(|2~T]) yields 

1 



G(T, R) = ( exp 



dt y z 



The integral is Gaufiian and can be done analytically 
given by 



VTy F 

VT yi 
For x\ = Xp 



G{T,R = 0) 



uoT 



(86) 

0, it is simply 

(87) 



sinh cuT 

In Fig. |Hl we plot the result of our Monte-Carlo calculation with iV = 2 14 ppl for uj = 1. 

The agreement with the analytical result is excellent, coinciding for more than 10 orders 
of magnitude with a very small error. However, the Monte-Carlo calculation is not free 
from problems, as becomes visible, e.g., for the ground state energy. The latter can be 
calculated from the standard identity 



lim 



T 



In Fig. we plot — log(G(T))/T versus T for large values of T, comparing the analytical 
result to the Monte-Carlo calculation with = 10 4 and 10 6 worldline ensembles. As is 
obvious from the figure, the Monte-Carlo approaches the analytical result for intermediate 
T, but ultimately departs from the analytical result from some larger value of T on. This 
point depends on the statistics of the Monte-Carlo. This phenomena is well known in the 
literature and is usually referred to as the "classical collapse" jHZj- 

We can cast the problem into our language by introducing the distribution of potentials, 

i.e., 

G(T, R = 0) = J dv 2 P{v 2 ) exp 1 
where the self-interaction potential is a local expression here: 



:U 2 T 2 V 2 



9) 



V 2 



dt y 2 . 



(90) 
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Figure 8: Propagator for the harmonic oscillator for x\ = xp = 0, for N = 2 14 ppl and 
uj = 1 and 10 4 worldline ensembles. The Monte-Carlo calculation (points) is shown in 
comparison to the analytical calculation of Eq. ()87j) (line). 



From Eq. ([89)1 . we observe that the limit T — > oo is controlled by v 2 —>■ 0. However, 
the distribution of potentials P{v 2 ) in this small v 2 limit is poorly determined, owing to 
the lack of statistics, P{v 2 ) — > 0. A typical distribution P(v 2 ) is displayed in Fig. ITU1 for 
N = 2 14 ppl. For any fixed number of worldline ensembles, there is always a minimum 
value of the potential, ^min, such that for any T larger than that determined by 



2rri2 

-UJ 1 f2min 



(91) 



the potential distribution P{v 2 ) is not well represented by the given statistical ensemble. 
This can be observed in Fig. where an increase of the size of the ensemble enlarges the 
region in T for which the Monte-Carlo calculation is valid. 

Our alternative method circumvents this problem by using an analytical parameteriza- 
tion of P{v 2 ) which has an appropriate behavior for low v 2 . For instance, the normalized 
trial function 



P(v 2 ) =A r < 



-(a/v2+bv2) 



AT' 1 = 2 A /^ 1 (2v / ^), 



(92) 



gives a reasonably good description of P{v2) and reproduces the large T behavior. Here 
a ~ 0.11 and b ~ 18.4 are constants determined by a fit to the potential distribution. 
From Eq. ()89|). we immediately get a numerical estimate for the propagator, 



G(T, R = 0) 



2b ir 1 (2 v /a6 + aw 2 T 2 /2) 



2b + lu 2 T 2 



Kx(2Vab) 



(93) 



9 The normalization of the kinetic term here is different from the main text, in order to make contact 
with the standard quantum mechanical normalizations. Both normalizations are connected by a simple 
rescaling of the worldlines. 
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Figure 9: — log G(T)/T for the harmonic oscillator. Shown is the exact analytical result 
for u — 1 (dashed line) and the Monte-Carlo results with 2 10 ppl for 10 4 and 10 6 worldline 
ensembles. 

such that the large T limit gives the ground state energy, 



which is in good agreement with the exact answer. 

D Analysis of the potential fit function 

Here we analyze some properties of the auxiliary function F( Pi )(x) introduced in Eq. (JSHJ) 
which are relevant for the mass renormalization as well as the weak- and strong-coupling 
behavior of the effective action for the scalar photon in Sect. 03 This function is defined 



with the pi being constant fit parameters and x = gT as in the main text. We observe that 
the exponent becomes large (but negative) for small as well as large v with an extremum in 
between. This suggest to evaluate the integral approximately by steepest descent in order 
to obtain analytic information about its behavior, e.g., for larger propertimes T. To be 
precise, we perform the approximation to the following order, 




(94) 



by 




(95) 




(96) 
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Figure 10: Potential distribution P(v%) for the harmonic oscillator for N = 2 10 ppl and 10 e 
worldline ensembles. Also shown is the fit to the proposed distribution. 



where 



f(y) = X V -\- — g [ln(f ) -p 3 ) 
^P2 



(97) 



such that the extremum ve which is a minimum of f(v) is determined by = /'(^e); the 
latter leads to a transcendental equation 



that can be brought to a standard form, 



p 3 - ln(v E ) 



vl 



pi e P3 x 



w e 



with w = ps — In ve- 



(91 



(99) 



This form, z = we w , is the defining equation of the "product logarithm" w(z) by which 
we can express the extremum of the integrand as a function of x: 



1 w(p 2 2 e P3 x) 

Ve{x) = -2 - 

pi 



X 



(100) 



One important property of the product logarithm is its asymptotic behavior for large argu- 
ment, w(z) — > In 2 for z — > oo, such that the extremum goes to zero for large propertimes 
x = gT. By contrast, the "inverse width" of the fluctuations around the extremum diverges 
for large propertimes and vanishing extremum, 



1 1 + p 3 - In ve 



vl 



(101) 
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such that the steepest-descent approximation is expected to be accurate in this limit. 
Putting it all together, the steepest-descent approximation of Eq. becomes upon in- 
sertion of Eqs. JSZJ, (MB and fpTTTfl into Eq. 



x exp < — ^ 
P2 



(gTy+O-lvl) 
1 ln(p 2 e p3 #T) 



In (3 ^ ; )-P3 



(102) 



Here we have already inserted the asymptotic form of the product logarithm for large 
argument. In this large- gT limit, the term ~ exp(— (l/p|) ln 2 gT) dominats the decay of 
F( p .}(gT), which is weaker than exponential but stronger than power-like. As a conse- 
quence, F^.^gT) does not contribute any finite part to the mass renormalization, due to 
the prescription specified in Eq. (|54|) . but guarantees the finiteness of the effective action 
Tqa[^4] in the massless-0 limit. 

Let us now study the opposite limit of small x = gT corresponding to a perturbative 
expansion. For this, the first exponential in the definition of F^(gT) in Eq. can be 
expanded and the v integral can be done order by order to give 

00 i(m+l)[(m+l)p§+2p 3 ] 



F {Pi) (gT) = V2^ Pm ~, i-9T) m . (103) 

m=0 

The mth order in this expansion corresponds to the I = m + 1-loop order in perturbation 
theory. For large m, the leading growth of the expansion coefficient is ~ exp(m 2 /2— mlnm) 
by use of Stirling's formula. Hence, this perturbative expansion gives an asymptotic series 
in the coupling. 

For comparison, we mention the analogous expansion for the simple- fit result, 

£^0 m T(m + i)r(« + 1) 

implying a weaker leading growth of the expansion coefficient ~ m a . This series is abso- 
lutely convergent for \gT / f3\ < 1. 
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